Coherence enhanced optical determination of the 229 Th isomeric transition 
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Abstract 

The impact of coherent light propagation on the excitation and fluorescence of thorium nuclei in a crystal 
lattice environment is investigated theoretically. We find that in the forward direction the fluorescence 
signal exhibits characteristic intensity modulations dominated by an orders of magnitude faster, sped-up 
initial decay signal. This feature can be exploited for the optical determination of the isomeric transition 
energy. In order to obtain a unmistakable signature of the isomeric nuclear fluorescence, we put forward a 
novel scheme for the direct measurement of the transition energy via electromagnetically modified nuclear 
forward scattering in a three-level A configuration. 
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While atomic clocks based on the microwave ground state hyperfine transition of idd Cs UJ] 
have revolutionized the field of time-frequency metrology, general agreement exists that the next 
two orders of magnitude in precision pose serious challenges [2J]. To circumvent this problem, an 
exciting alternative is to develop clocks based on a nuclear transition. The narrow line-widths of 
nuclear transitions and the isolation from external perturbations promise amazing stability. The 
optical transition from the ground state to the isomeric, i.e., long lived, first excited state of 229 Th 
was originally proposed by Peik and Tamm in 2003 y|] as a suitable candidate for a nuclear clock. 
Out of the entire known nuclear chart, the 229 Th ground state doublet consisting of the ground 
state 229 #Th and the isomeric state 229m Th offers a unique transition well within vacuum ultraviolet 
(VUV) with energy below 10 eV [4]. This transition is spotlighted also due to its tempting potential 
for testing the temporal variation of the fine structure constant |5fl> building a nuclear laser in 
the optical range [g], or providing an exciting platform for nuclear quantum optics and coherent 
control iM-llOll of VUV photons. A key obstacle on the way is the uncertainty about the transition 
energy: measurements up to date have been only indirect, suggesting an energy of ~ 7.8 ± 0.6 eV 

HO- 

One intensely pursued approach for a more precise measurement of the 229 Th isomeric transi- 
tion and the realization of a nuclear frequency standard involves nuclear spectroscopy of a macro- 



scopic number of thorium nuclei doped in VUV-transparent crystals 111 , 
ing density of up to 10 18 Th/cm 3 , crystals such as LiCaAlF 6 111 ill or CaF 2 [ 



1211 . Due to the high dop- 



120 offer the possibility 



to increase the nuclear excitation probability, an important requirement especially considering the 
very narrow radiative width of the transition, at present estimated at approximately 0.1 mHz 
Both LiCaAlF 6 and CaF 2 have large band gaps and present good transparency at the probable 
transition wavelength, such that the interplay with electronic shells in processes such as the elec- 
tronic bridge II 1311 . internal conversion 111411 or nuclear excitation via electron capture or transition 
fl 1 311 can be neglected. The disadvantage of the crystal lattice environment is that inhomogeneous 
broadening due to temperature-dependent shifts of the transition energy and inparticular spin- 
spin relaxation compromise the Rabi Jlfjl . Ramsey Jl^i 17 1 or hyper-Ramsey |18] interrogation 
schemes, commonly used in high-performance frequency standards. Instead, clock interrogation 
by fluorescence spectroscopy was proposed II 1211 . Significant suppression of the inhomogeneous 
broadening is expected as long as all nuclei experience the same crystal lattice environment and 



19Q. However, 



are confined to the Lamb-Dicke regime, i.e., the recoilless transitions regime IlllL 
together with the broad-band nature of the excitation, these very conditions lead to coherent light 
propagation through the sample and enhanced transient fluorescence in the forward direction, with 
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a speed up of the initial decay depending primarily on the sample optical thickness. These effects 



are well known from nuclear forward scattering of synchrotron radiation (SR) 



210 driving 



Mossbauer nuclear transitions in the x-ray regime, and from the interaction of atomic systems 



2411 . but have so far never been addressed in the context of the 



with visible and infrared light 1122 - 
229 Th isomeric transition. 

In this work we investigate the effect of coherent light propagation on the excitation and flu- 
orescence signal of the VUV isomeric transition in 229 Th. We find that the line-width of the 
light scattered in the forward direction is dynamically broadened by several orders of magnitude 
(~ 10 6 ) due to the coherent decay of the large number of 229 Th nuclei in the crystal and the strong 
dispersion close to the resonance frequency. While the dynamical fluorescence in the forward 
direction would therefore not be suitable for clock interrogation, it offers an enhanced and faster 
signal, advantageous for both a high signal-to-noise ratio and a shorter detection time. In order to 
exploit the full potential of the dynamical broadening we put forward a quantum optics scheme 
based on quantum interference induced by two coherent fields in a three level A-configuration as 
a novel way to identify the isomeric transition energy. The proposed setup reminding of elec- 
tromagnetically induced transparency (EIT) provides a clear signature for the excitation of the 
nuclear transition and enhanced precision in the optical determination of the transition frequency 
compared to a direct fluorescence experiment using only one field. 

The reduced transition probability of the nuclear magnetic dipole (Ml) 229 Th isomeric tran- 
sition has been evaluated theoretically to B(M1) ~ 0.032 Weisskopf units [6, 125 ] . This value 
corresponds to the very narrow radiative transition width of approximately 0.1 mHz. Since such 
narrow-bandwidth sources are not available in the VUV region, broadband excitation, with ei- 
ther synchrotron or VUV laser light, is envisaged to be employed in order to resonantly drive the 
isomeric transition. As long as elastic, recoil-free scattering of the incident light occurs, the contri- 
butions of all nuclei are spatially in phase in the forward direction and interfere coherently [26] . To 
the forward signal contribute not only the resonant frequency components, which are absorbed and 
reemitted, but also the non-resonant components, which experience dispersion. Thus the time evo- 
lution of the forward scattering response does not follow a natural exponential decay as expected 
for fluorescence involving a single- scattering event, but exhibits pronounced intensity modulations 
characteristic for the coherent resonant pulse propagation Il26ll . This modulation is known in nu- 
clear condensed-matter physics under the name of dynamical beat and for a single-resonance and 



a short 6 (t) -like exciting pulse has the form [^IzZllzsl] E(t) oc f exp _r/2 Ji( v / 4fr) where 



E is the transmitted pulse envelope, r a dimensionless time parameter r = t/t with t denoting 
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the natural lifetime of the nuclear isomer, £ the optical thickness and J\ the Bessel function of the 
first kind. The optical thickness is defined as £ = NaL/A 112 811 . where N is the number density of 
229 Th nuclei, a the resonance cross section and L the sample thickness. For a dopant density of 
10 18 Th/cm 3 and a sample with L = 1 cm, we obtain £ ~ 10 6 . 

From the asymptotic behavior of the Bessel function of first kind we see that for early times 
r <C 1/(1 + £) immediately after the excitation pulse (chosen as t = r = 0) the response field 
has the form E{t) oc £ exp~( 1+ ^ T / 2 , showing the speed-up of the initial decay. A parallel has been 
established between this faster decay and the superradiant decay in single-photon superradiance 
Ii29ll . For the case of the 1 cm crystal doped with 229 Th considered above, the speed-up would cor- 
respond to an enhancement of the decay rate by six orders of magnitude. This value significantly 
reduces the gap between the incoherent relaxation time of the excited state population (the natural 
radiative decay) and the short coherence time due to the crystal lattice effects, bringing inhomo- 
geneous and dynamical broadening on the same order of magnitude. At later times, the decay 
becomes subradiant, i.e., with a slower rate comparable to the incoherent natural decay rate due to 
destructive interference between radiation emitted by nuclei located at different depths in the sam- 



ple. A symmetric excitation of the sample would lead to the creation of a timed Dicke state [|30j] 
and a prolonged superradiant decay. So far this has been achieved for 57 Fe nuclei embedded in the 
center of a planar waveguide in the first measurement of the collective Lamb shift in nuclei [30]. 
One should however keep in mind that the dynamical beat behavior is by no means restricted to 
nuclear systems and the x-ray radiation wavelength, having been observed over a large wavelength 
range A ~ 1 A— 1 m in a variety of systems (see [26] and references therein). 

In the case of 229 Th doped in VUV transparent crystals both conditions for (i) recoilless, co- 
herent excitation and decay and (ii) broadband excitation are fulfilled for scattering in the forward 
direction. The incident laser or SR pulse duration is much shorter than the nuclear lifetime and 
provides broadband excitation. Due to the crystal transparency at the nuclear transition frequency, 
the main limiting factor for coherent pulse propagation, namely electronic photoabsorption, is not 
present in this case. We are interested in the initial superradiant response which occurs on a time 
scale much faster than the incoherent radiative decay and can facilitate the optical determination 
of the isomeric transition energy. As a consequence of the hyperfine splitting affecting the 229 Th 
nuclei in the crystal lattice, an analytical result for the scattered field intensity is not available and 
the numerical solution of the Maxwell-Bloch equations 113 ill has to be evaluated instead. This treat- 
ment allows also to take into account the inhomogeneous broadening occurring due to spin-spin 
relaxation in the crystal sample. We use in the following the quadrupole structure with hyperfine 
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level energies given by E„ 



*(<?)< 



7 



oo ) <Yzz 



[3m 2 - I is ( g )(I is (g) + l)]/[4J is(ff )(2/ is ( g ) 



1)] 



Q2I1 . where Qi S ( g ) = 1.8 eb (3.15 eb), with e the electron charge, is the quadrupole moment 
of the isomeric (ground) nuclear state, = —(100 — 200) is the antishielding factor and 



(1-T 



00 JYzz 



-10 18 V/cm 2 electric field gradient 



of 9 Th with the electric quadrupole splitting 



33 



3411 . Fig.QJa) shows the energy scheme 



120 of the ground and excited 229 Th nuclear 



states of spins I g =5/2 and J e =3/2, respectively. Here we adopt the recently proposed 229 Th:CaF 2 



crystal Ill2ll under the thermal equilibrium condition [ot], i.e., all the nuclei are equally populating 
the two ground states of m g = ±5/2 at sub-kelvin temperatures. 

Plain fluorescence spectroscopy of the isomeric nuclear transition is plagued by a small signal- 
to-background ratio 113 511 . For both determination of the exact nuclear transition frequency and for 
the nuclear clock interrogation schemes, this is a major issue, especially due to the weak coupling 
of the nuclear transition to the radiation field. According to current experimental data [35J], the 
main source of background photons is the a-decay of the 229 Th ground state with the total counting 
rate on the order of MHz and An spatial distribution. A possible solution is offered by coherent 
light scattering, which ensures that a directional signal collecting in the forward direction benefits 
of a more favorable signal-to-background ratio. However, impurities and color centers (which 
unfortunately may occur during the experiment as irradiation effects) can cause residual absorption 
and may interfere with the 229 Th optical spectroscopy even in the forward direction. This is why a 
clear signature of the nuclear isomeric transition 229fi Th— >■ 229m Th is desirable. In the following we 
put forward a quantum optics scheme designed to provide such a nuclear fluorescence signature 
by employing two VUV fields in a A-type three-level system in a setup reminding of EIT ||3 
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FIG. 1: (a) Quadrupole level splitting for the Th ground and isomeric states in the Th:CaF2 crystal 
lattice environment, (b) A left circularly polarized probe field drives the |||) — > |||) transition, (c) The 
electromagnetically modified forward scattering setup. The red thick arrow denotes the strong coupling 
field while the blue thin arrow shows the weak probe field. The detunings of the probe and coupling fields 
are denoted by A p and A c , respectively. 
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A radiation pulse denoted in the following as probe driving the relevant nuclear transition shines 
perpendicular to the nuclear sample, as shown in Fig.QJb). We consider a left circularly polarized 



weak VUV probe field driving the m e 



rrir. 



-1 magnetic dipole transition, where m e and 



m g denote the projections of the excited and ground state nuclear spins on the quantization axis, 
respectively. Additionally, as shown in Fig.QIc), we apply a strong linearly polarized continuous 



wave (CW) denoted as coupling laser field driving the m e 



transition. A A-type scheme 



is formed due to the combination of initial nuclear population and the selected field polarizations. 
The purpose of applying the coupling field is to create a dressed state which splits the nuclear 
resonance into a doublet [36y via the Autler-Townes effect [37|], giving rise to electromagnetically 
modified time spectra. The resulting beating can be then used as a specific signature for the 
detection of the nuclear fluorescence photons. 

We study the coherent pulse propagation for both setups in Fig. [TJb) (assuming f2 c = 0) and 



Fig. Ec) b y nu merica 



equations 



301 



ly determining the dynamics of the density matrix p via the Maxwell-Bloch 



390: 



in 

-d t Clp 



H,p 



+ Ps, 



(1) 



with the interaction Hamiltonian given by 
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h 

H -- 



and the decoherence matrix 
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In the equations above pjj. = AjA* k for {j, k} e {1,2,3} are the density matrix elements of 



p for the nuclear wave function 



Ai| ||) + A 2 |||) + A 3 1||) with the nuclear hyperfine 



levels shown in Fig. EJa). Furthermore, (031,032) = (-\/2/3, — 2/y/TE) axe the corresponding 
Clebsch-Gordan coefficients. The matrix p s describes the decoherence due to spin relaxation 
(73i; 732, 721) = 27r x (251, 108, 30) Hz III2I1 and T = 0.1 mHz |6|] denotes the spontaneous decay. 
The parameter rj is defined as 77 = g, and we consider a target thickness of L — 1 cm. Further 
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notations are fi p ( c ) for the Rabi frequency of the probe (coupling) fields which is proportional to 
the electric field E of the probe (coupling) VUV radiation 113 ll 13911 and c the speed of light. 

Since the electric quadrupole splitting of ground states can be experimentally determined via 
standard nuclear magnetic resonance (NMR) techniques Ill2n . we may set the two detunings of 
the coupling and probe fields to be identical but unknown initially, i.e., A c = A p = A. We 
numerically solve Eq. (OQ) with £ = 10 6 and Vt c with a laser intensity of 2kW/cm 2 and scan the 
region < A < 10 10 r. The time spectra for the setup in Fig.QJb) are presented in Figs.[2ta). 
In the absence of the coupling field, the probe only interacts with a two-level nuclear system, and 
the corresponding time spectra are less sensitive to the detuning A. For A = 10 9 r and 10 10 r, the 
probe field will not create any nuclear excitation and just passes through the crystal. The 229 Th 
nuclei will start to coherently scatter probe photons when A < 10 8 r, and with smaller detunings 
the time spectra become identical for A < 10 7 r. Essentially we could approach the wanted energy 
of 229m Th with the precision of 10 7 r, i.e. about 1 kHz, by measuring the signal photons scattered 
in the forward direction. As a disadvantage, since the spectra do not present any specific feature 
to confirm the excitation of the isomer, background from other unwanted electronic processes 
may present an identical scattering pattern at probe laser frequencies far away from the nuclear 
resonance. 

We turn now to the two-field setup presented in Fig. [He), which is more detuning- sensitive 
and provides specific identification features in the scattering time spectra. In Fig. (2b) we present 




FIG. 2: The time spectra of the forward scattered signal (a) without coupling field, (b) with coupling field. 
For all spectra the detunings A p = A c = A take values between < A < 10 10 r. The yellow filled 
area below the gray dashed double-dotted line delimits the region < A < 10 5 r. In (a) the values for 
A = 10 9 r and 10 10 r are smaller than the displayed scale. 
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our results for the electromagnetically modified forward scattering spectra. With the influence 
of Q c the spectra are much more sensitive to the detuning, and probe photons start to interact 
with the 229 Th nuclei already at A < 10 10 r. The enhanced sensitivity of the electromagnetically 
modified setup in Fig. [flc) is due to the detuning- sensitive dispersion relation which leads to a 
unique time spectrum for each combination of input probe field and coupling field strength and 
detuning. A comparison between Figs.^a) and (2b) shows that the coupling field has introduced 
additional beatings in the spectrum which can serve as a clear signature of the nuclear excitation. 
Furthermore, a significant advantage of our scheme is that the shapes of the spectra are not identical 
until A < 10 5 r. A fit of the theoretical and experimental curves can therefore be employed to 
determine the nuclear transition frequency. By scanning the detunings A, i.e., by varying the 
known frequencies of probe and coupling simultaneously, several A-dependent forward scattering 
time spectra can be measured. A fit with the theoretical curves in Fig. (2b) allows the determination 
of the detuning value and thus of the 229 #Th— > 229m Th transition frequency down to a precision 
of 10 5 T, i.e., 10 Hz. Once the isomeric transition frequency has been identified, the fluorescence 
clock interrogation scheme using the incoherent response emitted at an angle to the excitation 
pulse direction can be employed for the nuclear frequency standard. 

A brief estimate of the experimental counting rate for the signal and background photons in the 
forward direction is performed for the case of a 229 Th:CaF 2 crystal with the size 3 mmx 3 mm x 
10 mm and 229 Th concentration of 10 18 cm -3 . Given the 7880 yr half life of the ground state 229 Th, 
and the assumed 0.3 background photon per a-decay Q350 , the total counting rate of background 
photons (in An solid angle) is then 0.75 MHz, which could be significantly suppressed down to 1.8 
Hz by registering only signals within 1° x 1° in the forward direction. The background generated 
by the VUV pulse itself in the forward direction can be eliminated with a chopper or employing 
the nuclear lighthouse effect |40j]. Thus, 20 Hz would suffice as counting rate of the scattered 
probe photons in the forward direction, corresponding to a signal-background ratio larger than 10. 
On the other hand, the spin relaxation makes the coherences decay according to p%i ~ e~ 731< since 
73i ^> T 112D. The coherent effects should therefore be observed within the time scale of ms 
allowing 100 ~ 500 shots per second. 

A CW laser source in the VUV region is presently available only within limitations. The 
KBe 2 B0 3 F 2 (KBBF) crystals 114 ill have been successful in generating narrow-band VUV radiation 
via harmonic generation owing to their wide transparency and large birefringence necessary for 
phase-matched frequency conversion processes in this frequency region 11421 14311 . In particular, a 



quasi-CW laser with a 10 MHz repetition rate and 20 ps pulse duration has been achieved |43|]. A 
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CW coupling VUV laser at around 160 nm wavelength could also be generated by sum frequency 
mixing in metal vapors or driving a KBBF crystal with TkSapphire laser Il44l.l45n. The circularly 
pola rized probe laser on the other hand may be generated via nonlinear sum-frequency mixing 



|46J], or a harmonic of a VUV frequency comb 11471. 14811 around the isomeric wavelength. In the 



absence of a CW coupling field for the electromagnetically modified forward scattering scheme, 
one can use a moderate external magnetic field (B = 1 G) to split the single transition in Fig.QJb) 
into a doublet via magnetic dipole interaction. The magnetic-field-dependent quantum beat arising 
offers a signature that the nuclear isomer has been excited and presents reasonable sensitivity to the 
laser detuning. Numerical results for the LiCAF crystal for which the calculation of the magnetic 
hyperfine splitting is straightforward show that the detuning sensitivity goes as far as A = 10 7 , 
i.e., one could determine the nuclear transition frequency within 1 kHz with this method. 

In conclusion, the dynamical broadening of the nuclear signal in the forward direction can be 
used as efficient tool to determine the isomeric transition frequency. A two-field setup reminding 
of EIT provides both a signature of the nuclear excitation and enhanced sensitivity to the field 
detuning from the resonance. This scheme demonstrates a novel way to solve the question of 
identifying the isomeric transition energy and lays the foundation for developing nuclear quantum 
optics in 229 Th. 

The authors gratefully acknowledge fruitful discussions with G. Kazakov and Y. Nomura. 
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